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Abstract 



^ | Weak field approximate solutions in the A — >• limit of a model of de Sitter 

i 1 gravity have been presented in the static and spherically symmetric case. Although 

(3jT). the model looks different from general relativity, among those solutions, there still 

exist the weak Schwarzschild fields with the smooth connection to regular internal 
solutions obeying the Newtonian gravitational law. The existence of such solutions 
would determine the value of the coupling constant, which is different from that of 
the previous literature. Moreover, there also exist solutions that could deduce the 
' galactic rotation curves without invoking dark matter. 

in ' 
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1 Introduction 



In the 1970s a model of de Sitter (dS) gravity had been proposed [1-4]. In this model 
the Einstein-Hilbert action with a cosmological term could be deduced from a gauge-like 
action besides two quadratic terms of the curvature and torsion. The astronomical obser- 
vation [5, 6] on the asymptotically dS behavior of our universe has increased interest in 
the model as it may offer a way to deal with the dark energy problem [7]. If the Einstein- 
Hilbert term is required to be the main part of the gauge-like action, the cosmological 
constant should be large. The large cosmological constant may be canceled out by the 
vacuum energy density, leaving a small cosmological constant [4]. But it is difficult to 
explain why the large cosmological constant and the vacuum energy density are so close, 
but not exactly equal, to each other. On the other hand, if the cosmological constant is 
required to be small [1-3, 8], the quadratic curvature term would become the main part 
of the action. Note that in one of the field equations, the quadratic curvature term only 
contributes to the symmetric trace- free part. It is worth checking carefully whether the 
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model under this case could explain the experimental observations. Actually it has been 
shown [9, 10] that the model with a small cosmological constant may explain the acceler- 
ating expansion of the universe and supply a natural transit from decelerating expansion 
to accelerating expansion without the help of introducing matter fields in addition to dust. 
It has also been shown [11, 12] that all torsion-free vacuum solutions of this model are the 
vacuum solutions of Einstein's field equation with the same cosmological constant, and 
vise versa. Therefore, one may expect that the model with a small cosmological constant 
may pass all solar-system-scale experimental tests for general relativity (GR). 

However, it has been pointed out [9] that the energy-momentum-stress tensor of a 
spinless fluid in the torsion-free case of this model should be with a constant trace. Ques- 
tions would then appear such as, could the torsion-free Schwarzschild-dS (S-dS) solution 
be smoothly connected to internal solutions with nonzero torsion, or are there any S-dS 
solutions with nonzero torsion? In fact, the different dS spacetimes with nonzero torsion 
in this model have been obtained in [13, 14], but they are sill not the S-dS solutions. On 
the other hand, S-dS solutions with long-range spherically symmetric torsion have been 
given [15] in some special cases (not necessarily under the double duality ansatz [16]) of 
quadratic models of Poincare gauge theory of gravity, but our model does not fall into 
those special cases. 

We would like to firstly check the existence of the S-dS solutions with nonzero torsion 
in the weak field approximation. The Newtonian limit of general quadratic models in 
Poincare gauge theory of gravity has been calculated [17, 18] in the 1980's. In those 
calculations, quadratic terms in the field equations have been thrown away as usual. 
However, in the weak field approximation of our model, the quadratic curvature terms 
could not be easily thrown away, for the reason that they are the main parts of one of 
the field equations. In fact, we may let the cosmological constant be A — > 0, then only 
those quadratic curvature terms would appear in the limit of the field equation which 
contains the energy-momentum-stress tensor of the matter field. On the other hand, as 
those quadratic curvature terms are symmetric and trace-free, they would not appear in 
the trace part and the antisymmetric part of the field equations. The Newtonian limit of 
the trace equations has been recently analyzed [19], but a more complete analysis of all 
components of the field equations is needed. 

As was well known, Newton's theory of gravity meets great difficulties in the expla- 
nation of the flat rotation curves [20] of spiral galaxies. The most widely adopted way 
to resolve this problem is the dark matter hypothesis. But up to now, all of the possible 
candidates of dark matter (such as neutralino, axion, etc.) are either undetected or un- 
satisfactory. In the meanwhile, there also exist some models [21-24] which could deduce 
the galactic rotation curves without involving dark matter. We would like to explore the 
possibility of a new explanation for the galactic rotation curves from the dS gravity model. 

The paper is arranged as follows. We first briefly review the model of the dS gravity 
in section 2. In the third section, after dividing a field equation into its trace part, 
symmetric trace-free part and antisymmetric part, we attain the A — > limit of the 
model and calculate its weak field approximation in the static and spherically symmetric 
case. The weak field approximate solutions contain the weak Schwarzschild fields with 
nonzero torsion, which could be smoothly linked to regular internal solutions obeying 
the Newtonian gravitational law. The coupling constant is determined by the existence 
of such solutions. Moreover, solutions that could deduce the galactic rotation curves 
without invoking dark matter are also attained. Finally we end with some remarks in the 
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last section. 



2 A model of dS gravity 

A model of dS gravity has been constructed with a gauge-like action [1-4] 
S G = J L G = J K[-%x{F ah F ab )\ 



6 

^[RabcdR a C — p(R ~ ~p) + pS a b c S a C ] (1) 



in the units of h = c = 1, where k is a dimensionless coupling constant to be determined, 
and 

J r ab = {dA+-[A,A]) ab (2) 



or explicitly 



T A B ab = (dA A B ) a b + A A Ca AA c 



Bb 

2„a A „ 7-1 , 

(3) 



Rab a /3 — I e. a a A e^f, I S^ a 
—l~S/3ab 



is a dS algebra-valued 2-form and 



- ( J-C r o°" ) w 

is a dS algebra- valued 1-form. Here A,B... = 0, 1, 2, 3, 4 stand for matrix indices (internal 
indices) and the trace in Eq. (1) is taken for those indices. In addition, {e a a } is some local 
orthonormal frame field on the spacetime manifold and T a p a is the connection 1-form in 
this frame field, where a,b... stand for abstract indices [25, 26] and a, (3... = 0, 1,2,3 are 
concrete indices related to the frame field mentioned above. The curvature 2-form R a b a p 
and torsion 2-form S a a b are related to the connection 1-form as follows: 

R a b° 'p = (dT a /3) ab + T a ia A T 7 ^, (5) 

S a ab = (de a ) ab + TV A e\. (6) 



Moreover, 



p d p a _ d qc no c 

-^abc -Li-aba e c e /3 i " ab " ab^-a 



Rab — Racb C , R ~ 9^ R, 



ab ■ 



In fact, if spacetime is an umbilical submanifold of some (l+4)-dimensional ambient 
manifold and with positive normal curvature, then A a and Tab could be viewed [8] as 
the connection 1-form and curvature 2-form (in the dS-Lorentz frame) of the ambient 
manifold restricted to spacetime. Here, an umbilical submanifold means a submanifold 
with constant normal curvature, such as the dS spacetime which could be seen as an 
umbilical submanifold of a 5d Minkowski spacetime with positive normal curvature. A a 
could also be seen [27] as the Cartan connection of a Cartan geometry modeled on the 
dS spacetime and based on the spacetime manifold, with J 7 ^ the corresponding curvature 



3 



2-form. The Cartan geometry is a generalization of homogenous spaces with fibre bundle 
language, and one may refer to [27] for more details. In addition, it should be noted that 
3 /I 2 is identified [8] with a small cosmological constant A here, which is very different 
from the viewpoint of [4] where I is identified with the Planck length. The signature is 
chosen such that the metric coefficients are rj a p = diag(— 1, 1, 1, 1). 

The total action is 5* = Sm + Sq, where Sm is the action of the matter fields and the 
field equations can be given via the variational principle with respect to e a a , T a p a : 

8 2 

p(G a b + J^9ab) + \R\ 9ab — ^RacdeRb & + 72 1^1 9ab 

8 8 4 1 

—-pjS c daS c b + 7jV c S' a b c + -p^S ac dTb c H — S a fe = 0, (7) 
t z l z l z K 

4 8 1 

*->da 1 nr-pa r>de 



where 



Ta A V7 7~>da 1 r\rpa r>de C a , _ cl 

be — 4Vrfit be + ^-L de-ti be ~ J^^[bc] + ~Jbc 



Gab — Rab — -R9ab, T C ab — S° a b + 25 C [ a 5" i f ) ]rf, 



I p|2 p jjabed I C|2 c oabc 

l-fl-l — -Llabcd-tl j |>J| — "Jafc" , 

E a a = SSM/^e a a , S{, a = T> a a e a b, 

^- /3a XC /IF" ^- ca _ 8a„ce „ c 

ol = ob M /oL /3a, n = r a p e b ep 
and the variational derivatives are defined as follows: if 



5S M = J(X a a 5e a a + Y al3 a 5T^ a ), 



then 

5S M /5e a a = X a a , SS M /ST al3 a = Y[ aj3 f. 



3 Weak field approximation in the case with A — > 

As A is very small, it is interesting to see the case with A — > (/ — > oo). If A — > 
is directly set in the first field equation, then only the quadratic curvature terms are 
left, which are symmetric and trace-free. Thus, we would like to perform the following 
procedure. Divide the first field equation into its symmetric trace- free part, trace part, 
and antisymmetric part, then let I tend to infinity (I — > oo) in the above three parts and 
in the second field equation. The limiting equations are: 

\R\ 2 9ab — ^RacdeRb Cde = 0, (9) 

- R - V c S bc b + \s bcd T bcd + (/ 2 /8«:)S = 0, (10) 

R[ab] + V c S[ ab ] C + -S[ a Cd T b ]cd + (^ 2 /8«)£[a&] = 0, (11) 

- 2W d R da bc + T a de R de bc = 0. (12) 

When I — > oo, I 2 / k should tend to a finite value, otherwise Eqs. (10) and (11) would 
give X = and Euy = 0, which are unreasonable. In the torsion-free case, the scalar 



4 



curvature would be a constant from Eq. (12), and, therefore, £ = const from Eq. (10). 
This property has been pointed out by [9]. Now we are going to consider the week field 
approximation of the above equations. It would be assumed that 

9ab = Vab + lab, lab = O(s), S C ab = O(s), (13) 
£«6 = 0(3), T ab C = 0(3), (14) 

where s is a dimensionless parameter, called the weak field parameter. We will restrict 
ourselves to the static and 0(3)-symmetric case, with the static spherical coordinate 
system {t,r,9,ip}. r] ab could be defined by its components r]^ u = diag(— 1, 1, 1, 1) in the 
approximate inertial coordinate system {x^} is related to {t,r,0,<p} as usual: 

x° = t, x 1 = r sin 9 cos (p , x 2 = r sin 9 sin cp, x 3 = r cos 9. 

It could be proved [15, 25, 28] that i ab and S c ab have only these dependent components 
in the static spherical coordinate system: 

7oo = -20(r), 7 rr = ~2ip(r), (15) 

S° 0r = f(r), S r 0r = h(r), 

S e re = g(r), S e oe = -k(r), (16) 

S% = 0(r), S^ 0v = -k(r), 

where <fi plays the role of the Newtonian gravitational potential [1-3, 25, 29] and ip is 
an unknown function. It can be shown that components of i ab and S c ab in {x^} are as 
follows: 

7oo = -20, 7oi = 0, jij = (-2i))xiXj/r 2 , (17) 

S°0i = f x i/ r i = 0, 

S^j = (h + k^xj/r 2 - kStj, (18) 
S l jk = (-g/r)(5 l jX k - 5 l k Xj). 

For this case the contorsion tensor is related to the torsion tensor by 

K a bc = S cba . (19) 

Let T c ab = T a fJil/ d rT c (dx fl ) a (dx 1J ) b , where T a ^ is the connection coefficient in T c ab 
and the curvature tensor have the following first order approximate expressions: 

T c ab = ~ (dalb c + d bla c - d c lab ) - K c ab , (20) 

Rabc d = -(d c d [alb] d - d d d [alb]c ) + 2d [a K d m , (21) 

where 

K C ab = -j(S C ab + S ab c + S ba c ) (22) 

is the contorsion tensor. In the case with T> ab = p(r)(dt) a (dt) b , the equations for the first 
order approximate g ab and S c a b £1X6 clS follows: 

\R\ 2 Vab - AR acde R b cde = 0, (23) 
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R + d c S bc b + (l 2 /8n)p = 0, (24) 

R[ab] + d c S [ab] c = 0, (25) 

d d R da bc = 0. (26) 

The first order approximation of Eq. (9) is an identity. Here Eq. (23) is the second 
order approximation of Eq. (9). It should be considered for the reason that it is still an 
equation for the first order approximate g ab and S c ab . 

Now we are going to solve the above weak field equations. Applying Eq. (21) to Eq. 
(26), we have 

d d (d% clb] ^ + d^ d K cb ^) = 0. (27) 
The OOi (abc) component of Eq. (27) is 

d d (d d d aoo + 2d d K m ) = 0, 

which gives 

A[(<t>' + f) Xi /r]=0, 
i.e. (0' + f)"/r + 2(0' + f)'/r 2 - 2(0' + f)/r 3 = 0. 
Solving this equation we get that 

0' + / = Cr + D/r 2 . (28) 
The Oij component of Eq. (27) is an identity. The iOj component of Eq. (27) is 

d d d d K j0i - d%K j0d = 0, 

which gives 

6ij(-Ak - h'/r) + (xiXj/r 2 )[Ak + h'/r - 2k' /r - 2{h + k)/r 2 ] = 0, 

i.e. Ak + h'/r = 0, Ak + h'/r - 2k' /r - 2(h + k) /r 2 = 0, 
or equivalently 

h + k + rk' = 0. (29) 

The ijk component of Eq. (27) is 

d d d d d kl]l - d d d d d jlki - d d d l d kl]d + d%d jlkd + 2d d d d K kji - 2d d d i K kjd = 0, 
which gives 

A[(2Wr 2 - 2g/r)(S ikXj - S ijXk )} = 0, 
i.e. {2i\)/r 2 - 2g/r)" + (4/r)(2^/r 2 - 2g/r)' = 0. 
Solving this equation we get that 

tp/r 2 - g/r = B/r 3 + A. (30) 

From Eqs. (21), (17) and (18) the components of R a bcd i n {a^} could be attained as 
follows: 
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RoiOj = (0" + f')xiXj/r 2 + (0' + f)(5ijr 2 - x&^/r 
R ljk0 = (2/r 2 )(h + k + rk')x {l 6 j]k} 

Rijki = {^/r 2 - g/r)'(2/r)x[i(5j] k xi - 5j\ix k ) - 4(tp/r 2 - g/r)Si[ k Si]j. 

6 



Substituting Eqs. (28), (29) and (30) into the above equation, we get that 
Roioj = (Cr + D/r^'xiXj/r 2 + (Cr + D/r 2 )(5 ij r 2 - XiXj)/r 3 , 

Roijk — 0, 
RijkO — 0, 

Rijki = (B/r 3 + A)'(2/r)x [i (5 j]k x l - 5 j}l x k ) - A(B/r 3 + A)5 i[k S l}j . 



(32) 



Then 



\R\ 2 Voo - AR 0cde R cde = 12(C 2 - AA 2 ) + 2A(D 2 - B 2 )jr\ 
\R\ 2r 1oi — AR 0c d e Ri cde = 0, 

\R\ 2 Vij ~ 4=RicdeRj° de = (33) 



(AC 2 - IQA 2 - 16CD/r 3 - 32AB/r 3 + 16D 2 /r 6 - 16B 2 /r 6 )Si 
+ (A8CD/r 3 + 9QAB/r 3 - 2AD 2 /r 6 + 2AB 2 /r 6 )xiXj/r 2 



Roo — 3C, 

Roi = RiO — 0, 

Rij = ~{C + D/r 3 + AA + B/r 3 )5 tj + [3(5 + D)/r 3 ]xiXj/r 2 , 
R = _6(C + 2A), 



(34) 



Applying Eq. (33), the symmetric trace- free equation (23) could be solved, resulting in 
the following relations: 

C = ±2A, B = ±D, CD + 2AB = 0. (35) 

From Eq. (34) one could see that R[ a b] = and thus the antisymmetric equation (25) 
gives d c S[ a bf = 0. Components of d c S a b c are as follows: 



d c S 00 c = -f - 2f/r, 

d c Soi c = 0, , x 

d c S l0 c = x t [h'/r + 2(h + k)/r 2 }, {60) 
d c Sij c = [-2(g/r) -r(g/r)']5 ij + (g/r)'xiXj/r. 



Therefore, d c S[ a b] c = results in 

h'/r + 2(h + k)/r 2 = 0. (37) 

Combining Eqs. (29) and (37) we have 

(h-2k)' = 0, h + k = -r(h + k)'/3, 

h = 2k + C u h + k = C 2 /r 3 . (38) 
Now we turn to the trace equation (24). From Eqs. (36), (28) and (30) we have 

d c S bc b = -f'-2f/r + 2r(g/ry + 6(g/r) 

= 2ip'/r + 2ip/r 2 + A(p-3(2A + C). (39) 

Substituting Eqs. (34) and (39) into Eq. (24) results in 

- (A0 + 2ip'/r + 2ip/r 2 ) + 9(2A + C) = (l 2 /8n)p. (40) 
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Equation. (40) is an underdetermined equation for and if>. To solve it, more conditions 
are needed. Suppose that the external solution is the weak Schwarzschild field, i.e., 
= if> = —GM/r and the internal solution is regular and satisfies A0 = AirGp. Also, it 
is assumed that the internal solution could be linked to the external solution smoothly 
at r = Rs, resulting in a complete solution. For this complete solution, there would be 
C = -2A 



[Gm{r)/r 2 ]dr + 0(0) (41) 



with 



m(r) = / p = / 4irpr 2 dr, 

JB3(r) JO 

rR s 

0(0) = -GM/Rs - / [Gm(r)/r 2 ]dr, 
Jo 



and 



^ = [-Gm(r)/r] (1 2 /Q4ttGk + 1/2). 
Note that i/)(R s ) = -Gm(R s )/R s , thus, 

$ = -Gm(r)/r, (42) 

k = 1 2 /32ttG. (43) 

Equation. (42) is just the same as the corresponding case in GR. The result C = —2A is 
in accordance with Eq. (35). The torsion solutions will be given later in the more general 
case. 

Actually, Eq. (40) can also be solved with other supplementary conditions. For 
example, instead of assuming A0 = 4irGp, we may let 

A0 = 4vrG(p + p), (44) 

where p is an arbitrarily given function and does not contribute to the energy-momentum- 
stress tensors of matter fields. Fixing k by Eq. (43), then from Eqs. (35), (40) and (44) 
we have 

<f) = G {[m(r) + fh(r)}/r 2 }dr + 0(0), (45) 



with 



^ = --\ m {r) +m(r)/2] + -(2A + C)r 2 (46) 
r 2 



m(r) = / p = I 4irpr 2 dr, C = ±2A. 
Jm 3 (r) Jo 



'B 3 (r) 

As the internal solutions are regular, from Eqs. (28), (30) and (38), there should be 
B = D = 0, C 2 = 0. Therefore, the torsion solutions corresponding to Eqs. (45) and (46) 
are as follows: 

f = Cr- G[m(r) + m(r)]/r 2 , (47) 

g = (4A + 3C)r/2 - G[m{r) + m(r) /2]/r 2 , (48) 

h = -k= l -C l . (49) 
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The former case which is compatible with the Schwarzschild solutions corresponds to 
the special choice with C = —2A and p = 0. In fact, another choice of p could deduce the 
galactic rotation curves without invoking dark matter. To fit the galactic rotation curves, 
the dark matter density profile pdm has been given for many spiral galaxies, for example, 
see Refs. [30, 31], where the following choice is made: 

PDM = 2*G(r* + a*y (50) 

In Refs. [30, 31], the mass distribution of the galaxies is modeled as the sum of the 
bulge and disk stellar components and a halo of dark matter. The rotation curves are 
used to determine the two halo parameters a and a. In our model, we may just let the 
gravitational contribution p to take the same form as the dark matter density profile 
and suitably choose the integration constants, i.e., p = pdm an d C = —2A, without the 
inclusion of any real dark matter. For this case, Eqs. (45) and (46) is almost equivalent 
to the corresponding case in GR with dark matter. The parameters a and a can be 
determined in the same way as that in Refs. [30, 31] and, therefore, with the same values. 
For example, for the galaxy NGC 2841, a = 232 km/s and a = 11.6 kpc fit the rotation 
curve well, for the galaxy NGC 3031, o = 86 km/s and a = 2.0 kpc fit the rotation 
curve well, and so on. To say 'almost equivalent' but not 'equivalent', is because the term 
m/2 in Eq. (46) is different from rh, which should be the case in GR with dark matter. 
Fortunately, ip would not affect the rotation curves in the leading order of approximation, 
since the rotation velocity v c in a galaxy at a radius r in the approximation is given by 

v 2 (r) = r(d(p/dr). (51) 

For higher order approximations, more work is needed to be done and experiments with 
higher accuracy are needed to check the corresponding results. 

One may argue that p can not be specified a priori and could only be determined from 
observation. Actually, such kinds of functions also appear in other models which attempt 
to explain the galactic rotation curves without involving dark matter, such as Milgrom's 
modification of Newtonian dynamics (MOND) [21, 22] and the modified Newton's gravity 
in Finsler space [24]. What we can do now is to point out the geometrical meaning of p. 
From Eqs. (47) and (48), there is the relation: 

AnGp = ~^[r 2 (f - <?)]' - 3(44 + C), (52) 

which shows that the dark matter density could be directly related to the spacetime 
torsion. 

When I — > oo, I 2 / ' k should tend to a finite value, as was mentioned before. Obviously 
Eq. (43) satisfies this requirement. In [1-3, 8], the coupling constant is chosen to be 
— l 2 /64iiG, which is different from Eq. (43). For this case, there exists no regular internal 
solution which satisfies A0 = 4irGp and has a smooth junction to the Schwarzschild 
solutions. In fact, with this choice, the ratio of the coefficient of the Einstein term G a b 
to that of the matter term T, a b in Eq. (7) would be 1 : (— 87rG), just like the case in GR. 
But the role of the Einstein term in our model is different from that of GR. Actually, 
from Eq. (34), the Einstein term only contributes a constant term to Eq. (24), while the 
torsion term plays an important role in that equation. 
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Generally, the torsion tensor can be decomposed [32, 33] into three irreducible parts 
with respect to the Lorentz group: the tensor part, trace- vector part, and the axial vector 
part. For static and 0(3)-symmetric torsion, the axial vector part vanishes automatically, 
the tensor part satisfies / = 2g, h = 2k and the trace-vector part satisfies / = —g, 
h = —k. By Eqs. (47) and (48), if the torsion field only contains the tensor part, the 
matter density should be a constant: 

p = 3{2A + C)/2nG. (53) 

If the torsion field only contains the trace- vector part, then 

4p + 3p = 3(4A + 5C)/4ttG. (54) 

When p — 0, the matter density has to be a constant, too. 

4 Remarks 

The weak field approximation of the A — > limit of the dS gravity model has been 
calculated in the static and spherically symmetric case. The matter field is assumed to 
be a smoothly distributed dust sphere with finite radius. It comes out that if and only if 
k = 1 2 /32ttG, there exist regular internal solutions which satisfy A0 = AtrGp and have a 
smooth junction to the weak Schwarzschild fields. Recall that the main part of the action 
is the quadratic curvature term, not the Einstein-Hilbert term. The existence of the above 
solutions is of significance. The choice of the coupling constant here is different from that 
of [1-3, 8], where k = —1 2 /64itG. The choice in [1-3, 8] may due to a comparison between 
the dS gravity model and GR. But the Einstein term in the dS gravity model plays a 
different role from that of GR, as it only contributes to a constant term in the weak 
field approximate equations. Actually, the metric components and torsion components 
are closely related by Eqs. (28) and (30), such that the curvature tensor has to take a 
special form (32). Moreover, one may let C = A = and B = D = 0, then R a b c d = 0, 
i.e., the Weizenbock spacetime [34, 35] would be attained. On the other hand, the torsion 
tensor plays an important role in this model. In the weak field approximate solutions, if 
the torsion tensor only contains the tensor part, the matter density should be a constant; 
if it only contains the trace-vector part, Eq. (54) should be upheld. In particular, the 
matter density should be a constant in the torsion-free case. 

The trace equation (24) results in an underdetermined equation for the metric compo- 
nents and ip in the weak field approximation. Solutions with A0 = 4ttG(p + p) could be 
attained, which can explain the galactic rotation curves without the help of introducing 
dark matter. The geometrical meaning of p has been given by Eq. (52). It is a geometric 
quantity related to the spacetime torsion and can play the role of dark matter density 
though it is irrelevant to the energy-momentum-stress tensors of matter fields. To see 
the higher order behavior of the model, more work is needed to be done and experiments 
with higher accuracy are needed to check the corresponding results. 

Finally, it should be remarked that all the above results need to be reexamined in 
the case with A / 0. But we can conclude, at least, that if there exist regular internal 
solutions which are in accordance with the Newtonian gravitational law and could be 
smoothly extended to the weak S-dS fields, the coupling constant should be chosen as 
k = 1 2 /32itG. 
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